Irregular Magnetic Fields in Interstellar Clouds and Variations in the Observed 

Circular Polarization of Spectral Lines 



W. D. Watson 1 , D. S. Wiebe 1 ' 2 ' 3 , and R. M. Crutcher 2 
University of Illinois at Urbana- Champaign, Urbana, IL 61801 
w-watson@uiuc.edu 

o : 
o ■ 

<N ; ABSTRACT 

O . 

The strengths of magnetic fields in interstellar gas clouds are obtained through ob- 
servations of the circular polarization of spectral line radiation. Irregularities in this 
magnetic field may be present due to turbulence, waves or perhaps other causes, and 
may play an essential role in the structure and evolution of the gas clouds. To infer in- 
\jq • formation about these irregularities from the observational data, we develop statistical 

relationships between the rms values of the irregular component of the magnetic field 
and spatial variations in the circular polarization of the spectral line radiation. The 
irregularities are characterized in analogy with descriptions of turbulence — by a sum 
of Fourier waves having a power spectrum with a slope similar to that of Kolmogorov 
turbulence. For comparison, we also perform computations in which turbulent mag- 
netic and velocity fields from representative MHD simulations by others are utilized. 
Although the effects of the variations about the mean value of the magnetic field along 
the path of a ray tend to cancel, a significant residual effect in the polarization of the 
emergent radiation remains for typical values of the relevant parameters. A map of 
£> . observed spectra of the 21 cm line toward Orion A is analyzed and the results are com- 

pared with our calculations in order to infer the strength of the irregular component 
of the magnetic field. The rms of the irregular component is found to be comparable 
in magnitude to the mean magnetic field within the cloud. Hence, the turbulent and 
Alfven velocities should also be comparable. 

Subject headings: ISM: clouds — ISM: magnetic fields — polarization — stars: formation — 
turbulence 
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1. Introduction 

The breadths of spectral line radiation from interstellar gas clouds at radio wavelengths arc 
recognized to be much greater than would be caused by Doppler shifts resulting from thermal 
motions alone at the temperatures of the atoms and molecules in the gas. These breadths are 
believed to be an indication of disturbances (most likely a form of turbulence or waves) which 
are fundamental to the structure and evolution of interstellar gas clouds and star formation (e.g., 
Vazquez-Semadeni et al. 2000, and reviews in Lada & Kylafis 1999 and in Franco & Carraminana 
1999). Although the conclusions of earlier investigations that relate detailed information about 
the spectral lines specifically to the properties of turbulence have been questioned (see Dubinski, 
Narayan, & Phillips 1995), recent efforts along these lines are encouraging (Miesch & Scalo 1995; 
Lis et al. 1998). Of particular importance is the role of waves or turbulence in causing a magnetic 
pressure that acts parallel to the direction of the mean magnetic field lines in the clouds. This 
has long been viewed as a promising way to understand why the rate of collapse of interstellar 
clouds along the field lines is not more rapid than is indicated by the observational data. Magnetic 
disturbances have been thought to decay more slowly than ordinary disturbances, and hence more 
likely to persist and provide the component of the pressure that seems to be necessary (Arons 
& Max 1975). However, the premise that the magnetic disturbances are longer lived has been 
questioned in recent investigations (MacLow et al. 1998; Stone, Ostriker, & Gammie 1998; Padoan 
&; Nordlund 1999). In any case, irregularities in the magnetic fields are expected to be associated 
with the bulk motions which are reflected by the non-thermal breadths of spectral lines. At some 
level, these should be evident in the observational data. 

Information about the direction of the magnetic field in the interstellar gas has long been 
inferred from the linear polarization of starlight that results from selective absorption by dust 
grains. However, the polarizing power of the grains that absorb starlight apparently is greatly 
reduced in gas clouds of even modest densities (e.g., Goodman 1996). Though such polarization 
data are extensive, they do not thus seem to be sensitive indicators for the irregularities of the 
magnetic field in dense clouds. In contrast, there is no doubt that the magnetic field in interstellar 
clouds is evident in the circular polarization of spectral lines at radio frequencies through the 
Zeeman effect. Although such data are not so extensive as is desirable for extracting information 
on irregularities in the field, advances in the observational capabilities are continuing. It is thus 
useful to examine quantitatively the sensitivity of the circular polarization of spectral lines to 
plausible irregularities in the magnetic field in interstellar clouds. In particular, the polarization 
associated with a ray of radiation is a sum of contributions by the gas all along the path of the 
ray. The effects of deviations of the field from its mean value thus tend to cancel for the radiation 
that emerges from the cloud and to reduce its sensitivity as an indicator of irregularities. The 
goal of our investigation is to determine what information remains in the circular polarization of 
radiation that emerges from an irregular medium. Though not the emphasis of our investigation, 
an important challenge is to identify observational quantities for which the predictions distinguish 
between actual MHD phenomena and a simple generic description of the irregularities in the fields. 
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Finally, we note that the emission of linearly polarized radiation has been observed from grains 
in dense clouds. This linearly polarized radiation holds promise to become a valuable tool (e.g., 
Hildebrand 1996; Rao et al. 1998) as more sensitive and extensive maps are obtained. 

In Section 2, we describe how the irregularities are characterized and how the calculations 
of the circular polarization are performed. We focus on the quantity that commonly is obtained 
from the observational data — the single value B{ of the magnetic field that best fits the circular 
polarization (Stokes V intensity) when the Zeeman splitting is much smaller than the breadth of 
the spectral line. By comparing the statistical properties of the Bf which are inferred from the 
observations with the properties of the B{ which we calculate, we propose to gain information 
about the actual strength of the irregular component of the field within interstellar clouds. As an 
example, data for the 21 cm spectral line toward Orion A are analyzed in this manner in Section 3. 



2. Calculations 

The calculations here are restricted to the limit that ordinarily is relevant for interstellar gas 
clouds — that in which the Zeeman splitting is much smaller than the breadth of the spectral line. 
The optical depth t(v) for the intensity of the spectral line as a function of Doppler velocity v can 
then be expressed as, 

riy) = j f(s)exp(-(v-v t ) 2 /v 2 h )ds, (1) 

where the integral is along a straight-line path of length s through the medium in which the turbu- 
lent velocity v t varies along the path. The integral can be understood as a sum of the contributions 
from components of the gas along the line of sight where the bulk velocity of the component 
at location s is v t (s). The exponential arises from the Maxwellian form for the distribution of 
atomic/molecular velocities at any location and depends upon the thermal velocity w t h = \j2kT jra. 
All of the necessary information about the excitation, the abundance of the relevant species, and 
the atomic/molecular data is incorporated into the opacity function f(s). As measured by the 
Stokes V intensity, the net circular polarization is then related to the difference Ty{v) <C 1 between 
the optical depths associated with the two senses of circular polarization, 

t v (v) = -2pv- h 2 I f(s)B(s)(v - v t ) exp(-(« - v t ) 2 /v 2 h ) ds, (2) 

where p is a constant that depends upon the particular transition and B(s) is the component of 
the magnetic field that is parallel to the path (and hence to the line of sight of the observer) . Note 
that when B is constant along the path of the radiation, 

Ty = pBdr{v)/dv. (3) 

If the spectral line in the gas is being observed in absorption in front of a strong continuum source 
with intensity Iq as occurs toward Orion A in the analysis in Section 3, the observed intensity / is 
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given by 

/ = 7 exp(-r) (4) 

and Stokes- F is given by 

V = r v I, (5) 

so that if B were constant 

V = pBdl/dv (6) 

and B could readily be determined from the observation of / and V . Even though the magnetic 
field may not be constant along the path of the ray of radiation, a "best fit" (in the least squares 
sense) value Bf often is obtained from equation (3), 

B { = J (T V dT/dv)dv/ J p(dT/dv) 2 dv (7) 

or by the analogous fit to equation (6). 

If there are irregularities in the magnetic field, the best fit values for Bf will not be the same 
for spectral lines that are emitted from different locations across the surface of the cloud. Even 
at a single location, irregularities will usually prevent equation (3) or (6) from being satisfied by a 
single value for the magnetic field at all Doppler velocities within the observed spectral line. These 
variations in the inferred Bf are related, though only indirectly through the emergent Stokes-l^, 
to the variations in the magnetic field B(s) along the path of the radiation. The purpose of the 
investigation here is to provide a statistical relationship between the variations in the values of 
Bf that are inferred from the observational data and the amplitudes of the actual irregularities in 
the magnetic fields that are the cause of such variations. While more sophisticated and involved 
analytical tools might be imagined, we will limit our attention here to two basic quantities, (i) The 
standard deviation of the values of Bf that are inferred from the radiation emerging at the various 
different locations which comprise the surface of the cloud 

(Bf) sd = {(Bf) - (Bf) 2 }^ . (8) 

(ii) The average across the surface of the cloud of the residuals ((B 2 ^) 1 / 2 ) that remain when best 
fits are obtained for Bf at the various locations where at a single location 

B 2 es = J (t v -pBfdr/dvfdv/ J (pdr/dv) 2 dv. (9) 

Irregularities in the velocity and magnetic fields in the gas are characterized here in a conventional 
manner by the exponent in the power law distribution for the amplitudes of the Fourier components, 
by the root mean square (rms) of the velocity and of the magnetic fields, and by the smallest wave 
number (corresponding to the largest scale length) to which the power law distribution extends. For 
turbulence, this scale length represents the largest scale length at which disturbances are injected 
into the gas. The scale length implies a "correlation length" , which we also use in describing the 
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calculations. Representative configurations (or statistical "realizations" ) for the fields are then cre- 
ated by statistical sampling with Gaussian distributions for the Fourier amplitudes. Such methods 
are standard and are described, for example, in our own previous investigations (Wallin, Watson, & 
Wyld 1998, 1999) as well as elsewhere (e.g., Dubinski, Narayan, & Phillips 1995) where turbulent 
velocity fields are created. The turbulent magnetic fields are created in the same way (e.g., Wiebe 
& Watson 1998). For comparison purposes, we present a few examples of the realizations created by 
statistical sampling where the slope of the power spectrum is exactly that given by the Kolmogorov 
exponent. However, we focus our main attention on calculations in which the spectrum is somewhat 
steeper. This steeper spectrum leads to velocities and magnetic fields in coordinate space for which 
the slope of the correlation function is in better agreement with that calculated utilizing the fields 
from the available MHD simulations (see below). We have computed the "structure function" for 
the magnetic and velocity fields, as well as for the mass density, from these MHD computations. 
We find that the variation of the structure function (as defined by, e.g., Frisch 1995) corresponds 
best to (separation) 4 / 3 , whereas the Kolmogorov variation is (separation) 2 / 3 . The variation of the 
structure function is related to the variation of the power spectrum (e.g., Frisch 1995). We thus 
focus on a distribution for the power spectrum in Fourier space in our statistical sampling (fc~ 7 / 3 ) 
that is somewhat steeper than that of Kolmogorov (fc~ 5 / 3 ). We have verified that the structure 
functions which are computed using the fields obtained from the steeper power spectrum do agree 
well with those that are computed directly from the fields of the MHD simulations. We note that 
the (A; -7 / 3 ) dependence that we infer may not be a basic property of MHD turbulence, but may 
be a result of the numerical resolution in the MHD computations and the way in which the distur- 
bances are injected. The representative configurations for the velocities and magnetic fields that 
are created according to the foregoing statistical procedure thus provide a "generic" description for 
the irregular fields in terms of a minimal set of basic parameters. 

In detail, the foregoing generic description certainly is incomplete — most notably, it incorpo- 
rates no correlation between the velocities and magnetic fields and it contains no information on 
the variations in the density of matter. To obtain a quantitative measure of the influence of these 
deficiencies, we also perform computations in which these quantities are provided by the results of 
numerical simulations by others for time- independent, compressible MHD turbulence. We refer the 
reader to the papers by these authors for a detailed description of the calculational methods (see 
Stone, Ostriker, & Gammie 1998). Certain basic aspects of these calculations should, however, be 
summarized. The turbulence is driven by continuously injecting disturbances at a range of scale 
lengths centered on one-sixth the length of an edge of the computational cube. This determines the 
minimum wavenumber in the spectrum of the turbulence, and hence the correlation length for the 
velocities. In the way that we measure the correlation length (see below), the resulting correlation 
length corresponds to about one-twelfth the length of the edge of the computational cube. Three 
cases are considered according to the relative strength of the mean magnetic field -B avg and the 
thermal gas pressure P t h as described by the usual "plasma parameter" (3 = (87i\P t h/-B 2 vg ) = 0.02, 
0.2 and 2. For convenience, we label these as S, M, and W for strong, medium and weak, respec- 
tively. For all three cases, the rms turbulent velocity in each of the three orthogonal directions is 
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close to the value of 2.9c s in the weak (W) field case where the velocities are essentially isotropic. 
With increasing strength of the mean magnetic field, the gas becomes somewhat anisotropic so that 
the rms turbulent velocities parallel and perpendicular to the mean field are 3.4c s and 2.8c s , respec- 
tively, in the strong (S) case. Although we are limited to essentially a single ratio of the turbulent 
velocity to the sound velocity in the available MHD simulations, this ratio of approximately three 
is representative for relevant interstellar clouds. The chief difference between the three MHD cases 
is in the relative strengths of the mean magnetic fields which are in the ratio 1(W):3.2(M):10(S) — 
assuming that the same Pth is adopted in all three cases. The ratio of the standard deviation of 
the magnetic field in each of the three orthogonal directions to the mean magnetic field in these 
simulations is approximately 5/3 (W case), 2/3 (M case), and 0.15 [parallel] /0. 2 5 [perpendicular] 
(S case). These ratios can be viewed as depending mainly on the Alfvenic Mach number — the ratio 
of the three dimensional rms turbulent velocity to the Alfven velocity. The Alfvenic Mach number 
is 5(W), 1.6(M) and 0.5(S) in the simulations (the sonic Mach number is approximately five in all 
three simulations). The turbulent magnetic fields are then similar in the three the 
turbulent velocities. To relate the dimensionless results of the MHD computations to physical units, 
it is useful to recognize that there is freedom to choose the values for two additional scaling factors 
after (3 has been specified. These can be the sound velocity c s and the average matter density. 
Representative values for molecular clouds (e.g., Crutcher 1999) are (3 ~ 0.01 to 0.1 and include 3c s 
within the range of rms turbulent velocities that are observed. Their gas temperatures typically are 
10 to 30K, though the Orion A cloud that is analyzed in Section 3 is somewhat warmer (~ 100 K). 

Two idealizations will be considered for the opacity function f(s) — (i) f(s) will be taken as 
a constant, and (ii) f(s) will be taken as proportional to the matter density. The latter case is 
considered only in conjunction with the fields from the MHD simulations where we have information 
about the variations in the matter density. 

In Figure 1, we show the results for the standard deviations (-Bf) s d that are obtained from 
the computed Stokes-F as a function of the various parameters. For the statistically created 
fields, (-Bf)sd is independent of the average magnetic field in the gas. The results can thus be 
presented usefully in terms of the ratio (Bf) s ^/B ims where B rms is the rms of the irregular (or 
turbulent) component of the magnetic field in the gas. That is, B(s) = -B avg + i?i rreg ular for which 



the (128) 2 rays that emerge from the grid points on the surface and perpendicular to the surface 
of our computational cube consisting of (128) 3 grid points. Calculations also are presented when 
the intensity and Stokes- y for rays from separate areas of 8x8 surface grid points are first summed 
before the B{ are obtained from equation (7). The latter is intended to provide an indication of 
the influence of the finite angular resolution of a telescope. The calculations are performed as a 
function of the minimum wavenumber (our parameter k m \ n as defined in Wallin et al. 1998) at 
which the irregularities are injected into the medium. The correlation length is then computed 
from the structure function of the resulting fields in coordinate (that is, ordinary physical) space. 
We define the correlation length here as the separation at which the structure function "becomes 




Calculations are presented when Bf is determined separately for each of 
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flat". Due (presumably) to computational imprecision, the structure function does not become 
completely flat but has a wavelike variation with an amplitude of a few percent. We thus take 
the correlation length to be the separation at which the structure function becomes flat after this 
wave has been subtracted. In Figure 1, the quantity (i?f)sd/-E>rms is then presented as a function 
of the number iV corr of such correlation lengths along an edge of the computational cube. Other 
than those at a single correlation length corresponding to N corr = 12, all of the computations are 
performed with fields obtained from the power spectrum (fc~ 7 / 3 ) that is somewhat steeper than 
Kolmogorov. Within the range of turbulent velocities in our computations (see Figure 2) which 
covers the range for interstellar clouds, any dependence of (-Bf) s d/-Brms upon the rms value of the 
turbulent velocities is negligible. For iV corr < 12, the longest wavelengths in the turbulence are an 
appreciable fraction of the size of the cube and the purely statistical variations in (5f) S( j/-B rms from 
one realization of the cube to another (with otherwise identical parameters) are appreciable. This 
is indicated in Figure 1 by the broader lines. These lines are obtained by averaging (Bf) s ^/B rms for 
several statistical realizations. For iV corr > 12, the variations in (B^)^/ B rms from one statistical 
realization to another are negligible. 

In the right-hand panel in Figure 1, we present (-Bf) s d/ B rms obtained by utilizing fields from 
the MHD computations of Stone et al. (1998). These computations are available only for a single 
correlation length which corresponds to N COTT ~ 12 because the turbulence in these computations 
always is injected with the same distribution of scale lengths. In the strong field case, the correlation 
length is somewhat different in the directions perpendicular to the average magnetic field. For ease 
of comparison, the (-Bf) s d/-Brms based on the statistically created fields for N COTT = 12 in the left- 
hand panel are indicated by the two horizontal lines in the right-hand panel. In addition to the 
turbulent magnetic and velocity fields, there is a non-zero average magnetic field in the MHD 
simulations. This influences the medium, and potentially the results of our computations (a non- 
zero average magnetic field has no effect on the results of computations in the left-hand panel with 
the statistically created fields). As noted previously, observations suggest that molecular clouds 
fall between the medium and strong field cases in terms of the plasma parameter (3. There is no 
reliable way to incorporate variations in the mass density into the computations that utilize the 
statistical fields. Direct comparisons can thus only be made with the computations for the MHD 
fields where variations in the mass density are ignored (open symbols) in computing the optical 
depths. The computations in which f(s) in equation (1) is proportional to the mass density are 
indicated by filled symbols. For the computations with the finite (8x8) surface areas and when the 
optical depths do not depend upon mass density, (Bf) s d/B Tms ranges from approximately 0.25 to 
0.30 in the case of medium (M) field strength according to whether the line of sight is parallel or 
perpendicular to the direction of the average magnetic field. For the analogous correlation length 
(-^corr = 12), the computation with the statistically created fields is midway between these values. 
When the optical depths are assumed to be proportional to mass density, (Bf) s d/ B ims can be seen 
to range from 0.33 to 0.39 when the same fields from the MHD computations are utilized. 

In Figure 2, we show the ratio {(B 2 ^) 1 / 2 ) / B xms calculated with fields that are similar to 
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those used for Figure 1. Unlike the results in Figure 1, ((-B 2 ^) 1 / 2 ) / B rms does depend upon the 
magnitude of the turbulent velocities. Thus, the curves in the left-hand panel in Figure 2 are 
labeled according to the ratio of the rms turbulent velocity to the rms thermal velocity. In addition 
to the curves which are computed with the steeper power spectrum, calculations are again presented 
for iVcorr = 12 that utilize the Kolmogorov value for the exponent in the power spectrum (circles 
and crosses). As in Figure 1, the right-hand panel shows ((B 2 ^) 1 / 2 ) / B Tms computed with the MHD 
fields. The turbulent velocities of the MHD fields are slightly different from our choices for the 
curves in the left-hand panel. The results of the computations with the statistically created fields 
which are presented for purposes of comparison in the right-hand panel (horizontal lines) are thus 
computed with an rms turbulent velocity that is not exactly the same as that of any of the curves 
in the left-hand panel. The agreement between the results of the computations with the statistical 
and with the MHD fields is generally good. The calculations that are intended to reflect a finite 
(8x8) angular resolution coincide for the statistical and MHD fields when the optical depths do not 
involve mass density. When the optical depths depend upon mass density (specifically, when / is 
proportional to density), ((B 2 ^) 1 ^ 2 ) / B rms is seen to be increased for the 8x8 resolution from about 
0.24 to about 0.30-0.33 for medium fields. The increase is greater for the strong fields. Because 
B rms is much less than B avg in strong field case, we suspect that these MHD fields are not so 
representative of interstellar clouds where these fields probably are more similar in strength. 



3. An Application to Observational Data 

A statistical analysis such as the foregoing requires measurements of Stokes- F in the radiation 
from a large number of locations on the surface of an interstellar gas cloud. There is only one 
existing data set that meets this requirement — an H I absorption-line Stokes / and V map across 
the surface of a region toward Orion A (T. H. Troland, R. M. Crutcher, D. A. Roberts, W. M. 
Goss, k, C. Brogan, in preparation). The observations were performed with the VLA in the "C" 
configuration resulting in a spatial resolution of approximately 0.04 pc (15 arcseconds). The VLA 
synthesis mapping technique provides uniform sensitivity to Stokes / and V over the mapped 
region with no missing spatial points, which is a desirable characteristic for our analysis. Although 
observational data are available for a larger region, we limit the analysis to the region in which 
the background source is strong enough that the detections of the magnetic field are generally at a 
"three sigma" confidence level or better. With the angular resolution of the instrument, the region 
that meets our criteria for analysis consists of approximately 100 non-overlapping beams. Hence, 
we can obtain that number of independent observational values for Bf and B 2 es as a statistical 
sample for analysis. Within this area, there are no positions where the magnetic field was not 
detected at our sensitivity cutoff. Thus, the data set has no bias for positions with stronger 
magnetic fields or against positions with weaker magnetic fields. However, this data set does 
have significant limitations for an astrophysically meaningful analysis. Toward Orion A the 21 cm 
spectrum consists of a pair of partially overlapping absorption lines as indicated in the spectral 
data displayed in Figure 3 from a representative location on the surface. These two H I velocity 
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components arise in a feature which has been called the Orion lid. An extensive analysis of H I 
absorption due to the Orion lid has been carried out by van der Werf & Goss (1989) based on earlier 
Stokes I mapping. It consists primarily of neutral, atomic gas that probably is located close to 
and on our side of the H II region Orion A. The lid extends over at least 1.6 pc in the plane of the 
sky. Van der Werf & Goss argue that the more negative velocity H I is primarily photodissociated 
H2, some of which has been shocked, and that the more positive velocity component arises in the 
near envelope of the molecular cloud that is behind the H II region. We only utilize data from the 
positive velocity side of the intensity minimum that occurs in Figure 3 at approximately 5 km s _1 . 
This is primarily because the line shape of the other component leads to typically a factor of two 
lower sensitivity to Bf , so fewer positions are available for statistical analysis. In addition, the more 
positive velocity component is likely to be more relevant to the general interstellar medium, and 
we wish to avoid confusion in the analysis due to the overlap of these two components. Specifically, 
only spectral channels on the positive velocity side and where the absorption is less than ninety 
percent of the maximum absorption are included. The data give only the optical depth of the H I 
line. In order to infer column densities, one must know the spin temperature of the line. Toward 
the main exciting star of the H II region, 6 1 C Orioni, the H I column density has been measured 
with the ultraviolet Lyman-a absorption line; when compared with the 21-cm line optical depth, 
this yields a spin temperature of about 100 K. However, there is no information about possible 
variations in spin temperature over the mapped area. There is only the upper limit H2/H I < 10~ 4 
available; at least toward 6 1 C Orioni, the Orion lid is not a molecular cloud. Also, there is no direct 
information about the volume density of H I. An assumption of a spherical cloud and the fact that 
absorption is seen over at least 1.6 pc in the plane of the sky leads to «h ~ 10 3 cm~ 3 . However, 
it seems likely that the Orion lid has a more sheet-like geometry, which would raise the volume 
density, but by an unknown amount. Because of these limitations, our analysis of this Orion A data 
set should be considered to be primarily an example of the application of our analysis technique 
rather than one from which firm astrophysical conclusions should be drawn. 

The formal mean square error associated with the least squares fit to obtain Bf in equation 
(7) from a single, observed spectral line profile is 



where N is the number of intensities within the spectral line that are used to find the value of Bf 
from that spectral line. For the spectra included in Figure 3, the confidence level thus corresponds 
to Bf > 3a { . 

Differences between Stokes- V and dl/dv that lead to o\ can result from noise in the background 
source and in the instruments, as well as from the irregularities in the magnetic field that we are 
seeking to understand. In the relevant regime (ry <C 1 and V/I <C 1), the noise makes a contribution 
to d| that can be expressed as, 



a* = BlJ(N-l), 



(10) 




(11) 
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where Cy th is the mean square variation in the measured V that would occur even if the magnetic 
field were constant along the path of a ray. The sum is over the index i that designates the Stokes 
parameters Vi which are used in the computation of a particular value of Bf. The quantity 0y th 
is computed for an adjacent interval of velocities that is the same size as that used in computing 
Bf, but which is located outside of the spectral line. At these velocities, measurements of V are 
not influenced by variations in the magnetic field whereas the variations due to the noise should be 
essentially the same as within the spectral line. We thus utilize 

4th = E^ 2 /(iv-i), ( 12 ) 

3 

where the index j now designates locations in velocity that are spaced equivalently to those that 
are used in computing Bf, but which are located in the interval outside of the spectral line. The 
quantities of and of h are computed for each profile and are then averaged over the face of the 
cloud. We find for these data toward Orion A 

4>f * 0.1. (13) 

Since a^ h /af is relatively small, for simplicity we will ignore the contribution of the noise in our 
further discussions and make the simplifying approximation that <jf and (-Bf) s d are due entirely to 
variations in the magnetic field. 

From these observational data toward Orion A, we find (in agreement with Troland et al.) 
for the mean value of Bf obtained by averaging over the approximately 100 independent telescope 
beams that are included in our analysis 

(B f ) ohs = -170fiG, (14) 
for the standard deviation of these values 

(Bf)^ 8 = 130^G, (15) 

and for the average of the residuals 

<(£rcs) 1/2 }° bs = 83^G. (16) 
A histogram of the approximately 100 values of I?f obs used to find (Bf) ohs is given in Figure 4. 

To apply our calculations in Figures 1 and 2 for the purpose of inferring the strength of the 
irregular component of the actual magnetic field in the Orion A cloud as parameterized by B rms , the 
size of the actual gas cloud measured in correlation lengths N cori enters. Of course, real interstellar 
gas clouds are unlikely to be statistically uniform entities as are our idealized cubes. We can 
nevertheless estimate a correlation length from the variation of the measured quantities across the 
face of the Orion A cloud. That is, the structure function is computed for the 21 cm line and a 
correlation length is inferred in the same way as was done for our theoretical velocities in obtaining 
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Figures 1 and 2. In this way, we infer that there are approximately ten correlation lengths across the 
face of the region that is included in our analysis. The size of each telescope beam then represents 
about one correlation length. Hence, the relationship between the beam size and the correlation 
length is known, but the relevant size of the cloud (which is needed to obtain N corr ) is uncertain. 
Observations indicate only that the face of the entire cloud is greater than about 1.6 pc or about 
forty correlation lengths. The dimension of the cloud along the line of sight is unknown, but most 
likely is considerably smaller. Computations are thus performed for the quantities in Figures 1 and 
2 by first summing the intensities within non-overlapping surface areas that are one correlation 
length on a side. These are then averaged over the surface of the cube. Previously, we described 
first summing over fixed areas of the surface of the cube to obtain the curves for "8x8 rays" in 
these Figures. Now the sum is over a number of rays which varies with iV corr . This seems to be 
the most useful way to relate our calculations to the available observational data for Orion A. The 
results are represented by the lines labeled "Obs" in Figures 1 and 2. The "Obs" curve in Figure 
1 is somewhat less sensitive to the uncertainty about N corr than are the other curves. In Figure 
2, the "Obs" curve is quite insensitive to N COTT and has a value that is close to 0.4 over the entire 
range in Figure 2. The actual B rms that is then implied by the value of 83 (iG for ((-Brcs) 1 ^ 2 ) 0138 i s 

B Tms ~ 210/zG, (17) 

which indicates that -B rms ~ (B{) obs . The observed value {Bf)°^ s = 130/xG, together with the 
"Obs" curve in Figure 1, seem to imply a similar, though somewhat larger, value for B rms . The 
ordinate for the "Obs" curve ranges from approximately 0.35 at iV corr = 10 to slightly less than 0.2 
at iV corr = 50. If N corr > 10 as seems most appropriate, the inferred random magnetic field is 

B rms > 370^0. (18) 

In view of the involved procedures that are used in these two assessments of -B rms , the apparent 
discrepancy is not surprising. If the optical depths increase with mass density, the MHD calculations 
in Figures 1 and 2 indicate that the sense will be to raise the calculated values in both Figures. 
The increase will be a larger fraction in Figure 1 than in Figure 2 — an adjustment that is in the 
correct direction to reduce the discrepancy between the values for B rras that are inferred in the two 
approaches. An adjustment of this nature will also decrease the value for B ims that is inferred from 
((£rcs) 1/2 )° bs to about 130 fiG if one simply scales linearly from the results in the right hand panel 
in Figure 2. Our best estimates thus yield the conclusion that B rms is similar in magnitude to the 
mean value of the magnetic field {Bf) obs for this gas toward Orion A. 

We cannot exclude quantitatively the possibility that large scale changes in the average mag- 
netic field make a significant contribution to the variations that are inferred for the magnetic field 
in the Orion A cloud. An indication that such changes are small in directions along the line of sight 
is evident in Figure 3. The mean values of the magnetic fields are essentially the same in the feature 
on the left-hand side (the feature that we are analyzing) and in the feature on right-hand side of 
the profile. Since these two features probably are separated along the line of sight, we can conclude 
that the magnetic field probably does not change significantly in this direction. Deviations between 
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nearby channels — which presumably are due to changes over much shorter distances — tend to be 
greater than the difference between the mean values in the two features. The Troland et al. (in 
preparation) map of Bf does show evidence of some systematic gradients across the map, though 
these variations are not simple and do not appear to be large enough to dominate. Note that 
{(B^) 1 / 2 ) / B rms is sensitive to changes in the average magnetic field along the line of sight, but 
not to changes across the face of the cloud. In contrast, (B{) s d/B vms is sensitive to changes that 
occur in directions across the face of the cloud, but not in the direction along the line of sight. 

4. Discussion 

A number of deficiencies can be imagined in the simplifying idealizations of our calculations. 
Nevertheless, the conclusion appears to be firm that plausible turbulent magnetic fields are likely 
to have a significant effect on the magnetic field strengths that are inferred from the observation 
of the circular polarization of spectral lines at radio wavelengths. That is, their effect is unlikely 
to "average out". Clearly the exact ratio depends upon specific considerations outlined in the 
foregoing. One-third is a representative estimate for the ratio (Bf) s a/B rms of the rms deviations 
of the observationally inferred field to the rms of the irregular (or turbulent) magnetic field that 
actually exists in the cloud. As expected, there is a tendency for deviations from the average 
to cancel when the contributions are summed along the paths of the rays of radiation. When 
the opacity function f(s) is a function of location s because it depends upon matter density or 
other considerations, (-Bf) s d/ B ims will naturally tend to increase. This increase can be viewed as 
a result of effectively reducing the number of elements along the path that contribute which, in 
turn increases the scatter simply due to statistics. In MHD simulations, there should be some 
correlation between the matter density and the strength of the magnetic field. This also tends to 
increase (-Bf) s d/-E>rms> though we suspect that the effect is small based on our examination of the 
correlation between density and field strength. The results for the strong field (S) case in Figure 1 
might seem to be in conflict with this conclusion. We suspect, however, that the relatively larger 
values of (-Bf) s d/i?rms in this case are due mainly to the "statistical" effect of greater contrasts due 
to the anisotropic variations in the density. 

In a larger sense, the idealizations here are an incomplete description for star forming clouds 
because gravity is absent in these idealizations. Gravitational attractions will tend to cause clus- 
tering of matter, and hence more extreme variations in f(s) than considered in our calculations. 
In contrast, the results when f(s) is constant in Figures 1 and 2 suggest that purely statistical 
sampling provides velocity and magnetic fields that may be adequate to understand the circular 
polarization of clouds where gravity is not a major effect and when the finite resolution of the 
telescope is considered. This is not an altogether happy conclusion since it potentially limits our 
prospects for identifying the physical phenomena for which the non-thermal line breadths and 
variations in the magnetic fields are indicators. 

Our interpretation that the irregular component of the magnetic field in the Orion A cloud 
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is similar in magnitude to the average magnetic field (B rms = |(£>j? bs )|) implies that the magnetic 
field is quite disordered in this cloud. That the bulk motions of the gas are able to deform the field 
to this degree tends to indicate that the turbulent kinetic energy is at least comparable with the 
energy in the magnetic field, and thus that the turbulent velocity is at least comparable with the 
Alfven velocity. In numerical computations of MHD turbulence (including the MHD fields that we 
are using) , the energy of the irregular component of the magnetic field tends to be similar to the 
kinetic energy of the turbulence (typically 30 to 60%; Vazquez- Semadeni et al. 2000). For Alfven 
waves, these two components of the energy are equal. From this, it also follows that the turbulent 
and Alfven velocities should be comparable when -B rms = |(-B° bs )|. Padoan & Nordlund(1999) 
have emphasized that MHD computations in which the turbulent velocities are super- Alfvenic lead 
to predicted features for molecular clouds that are different in important ways from computations 
with lower turbulent velocities. In view of the uncertainties in our present analysis for B rms , we are 
unable to say whether the turbulent velocities are somewhat larger or somewhat smaller than the 
Alfven velocity. Alternatively, the ratio of the turbulent and Alfven velocities could be obtained 
from the spectral linebreadth in Figure 3 and the average field strength in equation (14) if the mass 
density were known reliably. However, as discussed in Section 3, the mass density is poorly known. 
A mass density of about 10 4 H-atoms cm -3 is required to yield an Alfven velocity that is equal to 
the turbulent velocity. As discussed in Section 3, the Orion A cloud probably is "sheet-like" . At 
least along one line of sight, its molecular fraction is small. It may not, thus, be representative of 
star- forming, molecular clouds. We note that in other molecular clouds where the matter density 
is better known, the velocities of the internal motions have been interpreted as equal to the Alfven 
velocities (Crutcher 1999). 

We are especially grateful to J. M. Stone for kindly making available the results of the MHD 
simulations. We also acknowledge valuable contributions by A. B. Men'shchikov in developing the 
initial computational methods, and by H. W. Wyld in a number of helpful conversations. This 
work has been supported by NSF Grants AST98-20641 and AST99-88104. 
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Fig. 1. — (left-hand panel) The ratio (Bf) s d/ B rms versus the number N COTr of correlation lengths 
across the edge of the computational cube — based on fields created by statistical sampling. The 
uppermost (solid) line is computed when each of the (128) 2 rays is treated independently. The 
lowest (dashed) line is computed from the (16) 2 values of Bf that are obtained when the radiation 
from areas that are 8 grid points on a side is combined before Bf is determined. The open circle and 
the cross represent computations in which the Kolmogorov exponent is used for the power spectrum. 
As discussed in the text, a somewhat steeper power spectrum is utilized in the other computations 
for this panel. The segments that are indicated by broader lines are especially dependent upon 
the particular statistical realization of the fields. In this region, the curves represent the average 
from several realizations. The line labeled "Obs" is intended to incorporate the finite angular size 
of the beam for the observational data which are analyzed in Section 3. (right-hand panel) The 
ratio (Bf) s d/ B rms computed with the fields obtained from MHD simulations for the cases in which 
the average magnetic field is weak (W), medium (M), and strong (S). This ratio is presented for 
computations in which the radiation propagates parallel or perpendicular to the direction of the 
non-zero average magnetic field, both for individual rays (lxl) and for areas of 8 grid points on a side 
(8x8) as in the left-hand panel. The opacity function f(s) is considered to be either independent of 
mass density (open symbols) or to be proportional to the mass density (filled symbols). To facilitate 
comparisons, the two horizontal lines in this panel indicate the values of (Bf) s d/ B rms taken from 
the left-hand panel for individual rays and for areas of 8x8 grid points. 
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Fig. 2.— (left-hand panel) The ratio ((£ r 2 cs ) 1/2 ) / B 

rms versus the number iV corr of correlation lengths 
across the edge of the computational cube — based on fields created by statistical sampling. The 
curves are labeled by the value of the rms turbulent velocity in units of v t h- The upper three (solid) 
lines are for the (128) 2 rays treated independently. The lower three (dashed) lines are computed 
for the (16) 2 values of -B 2 es that are obtained when the radiation within areas that are 8 grid points 
on a side is combined before i? 2 cs is computed. Other aspects of the Figure have meanings similar 
to their meanings in Figure 1. The broader segments of the lines indicate that, as in Figure 1, the 
results are sensitive to the particular statistical realization of the fields in this region. The open 
circles and the crosses represent computations in which the Kolmogorov value for the exponent is 
used for the power spectrum. As in Figure 1, a somewhat steeper power spectrum is utilized in the 
other computations for this panel. The lowest circle and cross are computed with a rms turbulent 
velocity that is similar to that of the MHD calculations. The line labeled "Obs" is intended to 
incorporate the finite angular size of the beam for the observational data which are analyzed in 
Section 3. (right-hand panel) The ratio {(B^ cs ) 1 / 2 ) / B rms computed for fields obtained from the 
MHD simulations for the cases when the average magnetic field is weak (W), medium (M), and 
strong (S). Symbols have the same meaning as in Figure 1. To facilitate comparisons, the two 
horizontal lines in this panel again indicate the value of this ratio when statistically created fields 
are used, and when the computations are performed separately for each of the (128) 2 rays and 
when the radiation from areas of 8x8 grid points is first combined. These are computed with a rms 
turbulent velocity (2.05 in our units) which is similar to that of the MHD computations. 
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Fig. 3. — (upper panel) The observed intensity (Jy per telescope beam), Stokes V intensity (Jy per 
telescope beam) and the derivative dl/dv (Jy per km s _1 per telescope beam) of the intensity from 
a representative location on the surface of Orion A that is included in our analysis, as function 
of Doppler velocity, (lower panel) The values for the magnetic field (//G) inferred from B = 
V/(pdI/dv) from the spectrum in the upper panel at several Doppler velocities. Error bars that 
correspond to one standard deviation are shown. 
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Fig. 4. — Histogram of the observationally determined values of i?f obs in the region toward Orion 
A that is being utilized for our statistical analysis. 



